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ABSTRACT 



We present the detection of carbon monoxide molecules (CO) at z = 1.6408 towards the quasar SDSS J160457. 50+220300. 5 using 
the Very Large Telescope Ultraviolet and Visual Echelle Spectrograph. CO absorption is detected in at least two components in the 
first six A-X bands and one d-X(5-0) inter-band system. This is the second detection of this kind along a quasar line of sight. The 
CO absorption profiles are well modelled assuming a rotational excitation of CO in the range 6 < T cx < 16 K, which is consistent 
with or higher than the temperature of the Cosmic Microwave Background Radiation at this redshift. We derive a total CO column 
density of N(CO) = 4 x 10 14 cirr 2 . The measured column densities of S i, Mgi, Znn, Fen and Sin indicate a dust depletion pattern 
typical of cold gas in the Galactic disc. The background quasar spectrum is significantly reddened (u-K ~ 4.5 mag) and presents 
a pronounced 2175 A dust absorption feature at the redshift of the CO absorber. Using a control sample of ~500 quasars we find 
the chance probability for this feature to be spurious is ~0.3%. We show that the spectral energy distribution (SED) of the quasar is 
well fitted with a QSO composite spectrum reddened with a Large Magellanic Cloud supershell extinction law at the redshift of the 
absorber. It is noticeable that this quasar is absent from the colour-selected SDSS quasar sample. This demonstrates our current view 
of the Universe may be biased against dusty sightlines. These direct observations of carbonaceous molecules and dust open up the 
possibility of studying physical conditions and chemistry of diffuse molecular gas in high redshift galaxies. 

Key words. cosmology:observations - galaxies: ISM - quasars: absorption lines - quasars: individual: SDSS J160457. 50+220300. 5 



1. Introduction 

Quasar absorption lines provide a powerful tool to detect and 
study gaseous baryonic matter at all redshifts in a luminosity- 
unbiased way. Large column densities of neutral gas are re- 
vealed by the damped Lyman-Q' absorption lines they imprint in 
the spectrum of background quasars. Because of the large neu- 
tral hydrogen column densities (N(Hi) <: 10 20 cm -2 ), similar 
to what is observed along Galactic lines of sight, and the pres- 
ence of metals at different levels of chemica l enrichment (e.g. 
iPettini et al.lll997tlProchaska & Wo"ifell2002l) , it is believed that 
a large fraction of Damped Lyman-a systems (DLAs) are lo- 
cated close to regions of star formation at high redshift (see, e.g. 
Wolf e et all I20051 for a review on the subject). 

The typical dust-to-gas ratio in DLAs, is generally less than 
one tenth of what is observed in the local interstellar medium 
(ISM) and only a small fraction (~10-15 %) of them show de- 
tectable amounts of mol ecular hydrogen dLedoux et all [2003; 
iNoterdaeme et al.. 2008a]). Even in these cases, the molecular 
fractions are small compared to what is seen in the Galactic ISM. 
It is therefo re likely that most of the DLAs probe only diffuse 
neutral gas (Petitje an et all 12000), In contrast, the detection of 
the cold (T ~ 10- 100 K), dusty and molecular gas, a fundamen- 
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* Based on observations carried out with the Ultraviolet and Visual 
Echelle Spectrograph, mounted on the ESO Very Large Telescope, un- 
der Prgm. ID. 081.A-0334(B). 



tal ingredient for star-formation, is elusive till now in absorption 
studies. The corresponding regions could have been missed so 
far due to the large extinction they are expected to produce. It 
could also be that molec ular gas escapes detection be cause of its 
very small cross section dZwaan & Prochaskal 12006). 

Gamma Ray Bursts (GRBs) are thought to be located in 
star forming regions, and the associated absorptions are there- 
fore more likely to arise from the densest part of the ISM at 
high redshift. However, the physical state of the absorbing gas 
is very likely to be influenced by the intense UV radiatio n field 
arising from the GRB itself (e.g. IVree swiik et all 120071) . This 
could explain the absenc e of H2 in the majority of GRB -DLAs 
dTumlinson et all 120071) . However, only a handful of GRB- 
DLAs have been studied so far and existing data are still con- 
sistent with the statistics of H2 detections in intervening QSO- 
DLA sam ples (Ledoux et al. 2009, submitted). When molecules 
are seen (Prochask a et all l2009h . their excitation is observed to 
be high, indicating strong UV pumping from the GRB afterglow. 

We have started a programme to search for cold gas along 
QSO lines of sight. The huge number of quasar spectra available 
in the Sloan Digital Sky Survey (SDSS) and the faint magni- 
tude limit achieved by the survey allows us to identify absorp- 
tion systems having unique characteristics. In particular the sys- 
tems can be selected on the basis of the presence of C 1 which 
should flag predominantly cold neutral gas. We searched for 
strong (W r (C lA 1656) ~ 0.5 A) C 1 absorbers at z ~ 1 .5 - 3 along 
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the line of sight of ~ 40000 QSOs from the SDSS Data Release 
7. 

The selection lead to the first dete ction of carbon monox- 
ide (CO) absorption lines at z ~ 2.4 dSrianand et all |2008b). 
Similarly, selection of strong Mg n systems at intermediate red- 
shift lead to the detection of the 2175 A feature and 2 1 -cm ab - 
sorption in two z ~ 1.3 Mgn systems dSrianand et all 12008a). 
Till now, direct signatures of dust at z > 0, such as the UV bump 
or diffuse interstellar bands (DIBs) have only been reported in a 
small number of cases (e . g..lMotta et al.L [2002h fWucknitz et al.L 
120031; lEllison et all 120081: Liang & Lill2009t) . 

In this paper we present the simultaneous detection of carbon 
monoxide absorption lines and 2175 A dust feature in an inter- 
vening absorber at z = 1.64 towards SDSS J 160457 +220300. 
We present high spectral resolution observations in Sect. 2, the 
metal content of the system in Sect. 3 and the analysis of molec- 
ular lines in Sect. 4. We discuss the extinction and presence of a 
UV bump in Sect. 5. Finally, we conclude in Sect. 6. 

2. Observations 

The quasar SDSS J160457.50+220300.5 (z em = 1.98) was ob- 
served with the Ultraviolet and Visual Echelle Spectrograph 
(UVES; iDekker et all 120001) in visitor mode on June 28, 29, 
and 30, 2008. The total exposure time on source is 29400 s. 
Both UVES spectrographic arms were used simultaneously tak- 
ing advantage of a dichroic setting with central wavelengths 
of 390 nm in the blue and 564 nm in the red. The result- 
ing wavelength coverage is 329-451 nm and 462-665 nm with 
a gap between 559 and 568 nm corresponding to the phys- 
ical gap between the two red CCDs. The CCD pixels were 
binned 2x2 and the slit width adjusted to 1" to match the 
ambient seeing conditions. This yielded a resolving power of 
R = 47 500 in the blue and R = 45 000 in the red as mea- 
sured from Th-Ar lines from the calibration lamp. The data 
were reduced using the MIDAS-based UVES pipeline v 2.9.7, 
which performs an accurate tracking of the object while sub- 
tracting the sky spectrum at the same time. Cosmic ray im- 
pacts and CCD defects were rejected iteratively. Wavelengths 
were rebinned to the vacuum-heliocentric rest frame and indi- 
vidual exposures were co-added using a sliding windows and 
weighting the signal by the signal-to-noise ratio in each pixel. 
We analysed the spectrum using standard Voigt-profile fitting 
techniques. Oscillator streng ths and wavelengths of C O absorp- 
tion lines were taken from Morton & Noreau (1994) with up- 
dated values from Eidelsberg & Rostaj d2003jT for the inter- 
band systems. Heavy eleme nt abundances are given relative to 
solar dGrevesse et all l2007h . with [X/H] = log N(X)/N(H) - 
log (X/H) and assuming N(H) = N(H i). For short reference 
through the paper, we refer to SDSS J160457. 50+220300.5 as 
Q 1604+2203. 

3. Metal content 

Absorption lines from singly ionised metals (Fe n, Si n, and Zn n) 
as well as neutral species (S i, Mg i, and C i) are detected around 
Zabs = 1-6405 over ~200 kms _1 (see Fig. [TJ. The covered ab- 
sorption lines of neutral species are exceptionally strong, with 
C 1/1/11560,1656 being heavily saturated and S i clearly detected 
in several transitio ns. Note that S i is rarely seen in QSO ab- 
sorbing systems dOuast et alll2008tlSrianand et al.Ll2008bl) . Fei 
is not detected down to log N(Fe i) = 1 1 .4 for each component 
(3 cr upper limit). We performed simultaneous multicomponent 



Voigt-profile fits with fitlyman dFontana & Ballesteri, 1 19951) to 
constrain redshifts, Doppler parameters b and column densities, 
see Table Q] Five narrow components are required to fit the neu- 
tral species (S i and Mgi). The profiles of singly ionised species 
(Zn ii, Si ii and Fe n) require additional broad components. 

While the narrow components can still be seen in the profiles 
of Zn ii transitions, the presence of additional broad components 
introduces a degeneracy in the results in particular the relative 
column densities in the narrow and broad components. We have 
somewhat artificially associated the broad and narrow compo- 
nents so that one could consider the sum of the column densi- 
ties in the broad (b > 10 kms -1 ) and narrow (b < 10 kms~') 
components to be representative of the column density in the 
six components listed in Table 1. Indeed, narrow components 
can be completely lost into the broad ones as it is the case for 
Fen and Sin. However, because lines are in the optically-thin 
regime, integrated column densities will not depend upon the 
actual number of components. The above decomposition shows 
that the coldest gas, as traced by the S i components, is found 
in clumps embedded in a more turbulent and probably warmer 
phase. This structure is natural (see e.g. iPetitiean et all 11992 ) 
and is usually observed in the H2 phase of DLAs dLedoux et al , 
12001) . 

The immediate consequence of this is that the depletion fac- 
tor estimated from integrated column densities may not be a cor- 
rect representation of the actual depletion factors in individual 
components. 

Using the integrated column densities (Table [TJ and upper 
limits of different ions (logA^(Nin) < 13.2 and logMCrii) < 
12.8 at 3 cr c.l., for the whole profile) we find, [Fe/Zn] = -1.47, 
[Si/Zn] = -1.07, [Ni/Zn] < -1.5 and [Cr/Zn] < -1.3. This is 
consistent with the depletion pattern seen in the cold neutral 
medium o f the Galacti c ISM and that of the Large Magellanic 
Cloud (see lWeltv et allfl999T) . 

From the equivalent width of C 1/11280, which is the weak- 
est Ci line available we derive logA^(Ci) > 15.3, in the linear 
regime limit. However, even this line is saturated, and filling the 
Ci profiles with five components with the same parameters - 
Doppler parameters and redshifts - as the S 1 components gives 
log N(C 1) > 16.0 which is still a lower-limit. Therefore, this sys- 
tem has by far the highest column densities of C 1 and S 1 known 
in any QSO absorption system. 



4. Carbon monoxide 

Carbon monoxide absorption lines are detected in several A-X 
bands and the d-X (5-0) inter-band system. 

Following [Burgh et"ail d2007l) . we use the excitation tem- 
perature as an external parameter. We fitted the CO p rofile us- 
ing a IDL code based on MPFIT dMarkwardA 120091) . which 
performs ^-minimisation by Levenberg-Marquardt technique. 
Components associated to the reddest S 1 components, at z = 
1.64071 and 1.64083, are clearly detected in the J = and 
J — 1 rotational levels (see Fig. [2]). The S/N ratio of the data 
is not high enough so that an independent fit can be performed. 
In addition, lines from different rotational levels are blended. We 
therefore fixed the redshifts and Doppler parameters to those ob- 
tained from S 1. We varied only the total CO column density for 
each of the two components, distributed among rotational levels 
up to J = 3, using a single excitation temperature T ex . 

If we assume the excitation of CO is due to the Cosmic 
Microwave Background Radiation alone, i.e. T ex = Tqmbr(z = 
1 .64), then the best fit model is achieved with = 1 .08 for a 
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Comp." 


z 


i(kms-') 


Znn 


Fen 


logJV(X) (cm- 2 ) 
Sin Mgi 


S i 


CO 


IB 


1.63967 


17.8 + 0.8 


12.31 ± 0.03 


13.73 ± 0.02 










2N 


1.63992 


6.3 ± 0.4 


12.34 ± 0.05 






12.85 ± 0.04 


13.18 ±0.02 


<13.64 


2B 


1.63996 


11.4 ±0.3 


12.09 ±0.11 


14.10 ±0.01 


14.72 ± 0.03 








3N 


1.64028 


2.8 ± 0.6 


12.07 ± 0.04 








12.64 ± 0.05 


<13.28 


3B 


1.64034 


37.6 + 1.4 




1 4 00 + 0? 


14.52 ± 0.09 








4N 


1.64080 


4.5 ± 0.4 


12.06 ± 0.05 


13.29 ± 0.04 


13.96 ±0.13 


12.73 ± 0.03 


13.01 ± 0.02 


<13.80 


5N 


1.64069 


6.8 ± 0.3 


12.70 ± 0.03 


14.08 ± 0.01 


14.45 ± 0.05 










1.64071 


3.6 ± 0.4 








12.66 ± 0.04 


13.10 ±0.02 


14.24 ± 0.20 


6N 


1.64083 


3.9 + 0.3 


12.42 ± 0.02 


13.67 ±0.02 


14.02 ± 0.10 


12.33 ± 0.07 


13.13 ±0.02 


14.34 ± 0.10 


total 






13.30 ± 0.02 


14.67 ±0.01 


15.13 ±0.03 


13.28 ± 0.02 


13.75 ± 0.01 


14.59 ± 0.11 b 



" Narrow components (b < lOkms ') are indicated by "N" while "B" stands for broad component (b > lOkms '). 

b Total N(CO) measured from comp. 5 and 6 only. Undetected components (2N, 3N and 4N) could increase this value by up to 0.1 dex (see text). 



total column density logA^(CO) = 14.6. The model is superim- 
posed on the observed spectrum in Fig. [2] 

J ex and 7cmbr are expe cted to be equal for low gas pres- 
sure (Sri anandetdl l2008b) and in the absence of UV pump- 
ing. Excitation by photon t rapping becomes significant only 
when N(CO)> 10 16 cm" 2 dBurgh et all 120071) . Unfortunately, 
unlike in the case of the z a b s = 2.42 CO-bearing system to- 
wards SDSS J143912+1 1 1740, all observed neutral carbon lines 
are saturated and it is therefore impossible to estimate the pres- 
sure of the gas from the population of C i fine-structure levels. 
We fitted different models along a grid of excitation temperatures 
and Doppler parameters. A good fit (y 2 S 1 .2) is achieved for a 
CO excitation temperature equals to or larger than the CMBR 
temperature and for b > 0.4 km s" 1 . We find 6 < T ex < 16 K at 
the 5 cr confidence level. 

Assuming a single excitation temperature for all rotational 
levels could be considered a rough assumption. However, Tqj is 
generally found co nstant for J < 3 in diffuse molecular gas in 
the local Universe (She ffer et all [2008). Moreover, most of the 
optical depth of the CO profile is due to absorptions from / = 
and J — 1 . Actually, a two rotational levels model gives the same 
results. We checked that the total CO column density only varies 
by less than 0. 10 dex regardless of the b and T ex values as long as 
the condition x 2 < 1-2 is fulfilled. The decomposition between 
the different components is quite uncertain however. 

CO could be also present in the three bluest S i components, 
although not detected directly at the 3 cr level. Because these 
components are far enough from the two detected CO compo- 
nents, they have very little influence on the excitation tempera- 
ture. However, they could contribute to increase the overall CO 
column density (over the range v = -150 to +50 kms~') by 
0.1 dex. 

From the lower limit on N(Ci) we derive N(CO)/N(C i) < 
0.05 which i s typical of what is me asured in the diffuse molecu- 
lar medium (Federm an et all[l980l) . 



5. Dust and 2175 A UV bump 

The presence of dust grains influences the physical state of the 
gas through photo-electric heating, UV shielding, and formation 
of molecules on the surface of grains. It can be deduced from the 
depletion pattern of different elements (see Sect. [3} and/or from 
the reddening of the background QSO. 

The spectral energy distribution (SED) of Q 1604+2203 is 
very red with u-K~4.5 mag. We note this value is only indicative 
as u (SDSS) and K (2MASS) magnitudes were measured with 4 
years interval and the exact u-K value could be different due to 
quasar optical variability. The flux-calibrated SDSS spectrum of 
Q 1604+2203 is visibly affected by reddening and also shows a 
clear curvature around 2175 A in the rest frame of the z = 1.64 
absorbing system (see Fig. [3). We investigate the possibility that 
the red colour of the quasar is due to the presence of dust in the 
CO-bearing system. 

We performed a ^-minimi sation between t h e data and 
a SDSS composite spectrum (IVanden Berk et all 12001 1) red- 
dened by different types of extinction curves, namely that from 
the Small Magellanic Cloud (SMC), Large Magellanic Cloud 
Supershell (L MC2), Large Magel lanic cloud (LMC) and the 
Galaxy (MW) dGordon et alll2003l) . The strength of the 2175 A 
bump is very different from one extinction curve to the other: 
while the UV bump is absent in the SMC extinction curve it is 
strongest in the MW and LMC extinction curves. In addition, the 
SMC extinction curve has a large UV extinction at A < 2000 A 
compared to the LMC and MW extinction curves. The average 
extinction curve of the LMC supershell (LMC2), which is part 
of the 30 Dor star forming region, has a bump strength and a UV 
extinction in-between the two extremes. 

Results of our best fit models using different extinction 
curves are presented in Tabl e [2] A detailed des cription of the 
procedure can be found in ISrianand et all d2008al) . The best re- 
duced x 2 value is reached for the LMC2 Supershell dust extinc- 
tion curve with E(B-V) = 0.27+0.02 for Ry = 2.7. 

We performed a test on a sample of SDSS spectra to rule out 
the possibility that a peculiar intrinsic spectral shape is mistaken 
for a 2175 A bump dPitman et all 120001) . Our control sample 
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Fig. 1. Voigt-profile fit to neutral and singly ionised species. 
The origin of the velocity scale is taken at z — 1.6408. Vertical 
dotted lines indicate the position of narrow components ('N' 
in Table [TJ. The contribution from the broad components only 
('B' in Tabled is represented by the blue profile. The extremely 
strong C i lines, though not fitted, are also shown in the top pan- 
els. The C 1/11280 profile has been smoothed by 2 pixels for pre- 
sentation purpose only. 

Table 2. Results of fitting the SED of Q 1604+2203 with differ- 
ent extinction curves 



Extinction law 


R v 


E(B-V) 


xt 


MW 


3.1 


0.26 


±0.01 


3.0 


LMC 


3.4 


0.24 


±0.02 


1.6 


LMC2 


2.7 


0.27 


±0.02 


1.2 


SMC 


2.7 


0.16 


±0.02 


2.1 



consists of all 529 SDSS quasars with emission redshift within 
Az=+0.01 of that of Q 1604+2203 (z em = 1.979). The distri- 
bution of S/N ratios of the spectra^ is shown in the upper right 
panel of Fig. [4] We fitted all quasar spectra using the SDSS com- 
posite spectrum reddened by a SMC extinction curve shifted 
to Zabs = 1.64, which does not present any UV bump. Using 
X 2 minimisation we estimated, for each quasar, E(B-V) and the 
strength of a possible 2175 A feature. The latter is calculated as 
the integrated flux difference between the fitted composite spec- 



1 Median signal-to-noise ratio in the i-band (SN_I) provided by SDSS 
together with the ID spectra (spSpec*.fit) 
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Fig. 2. Voigt profile fit to detected CO absorption lines. The ori- 
gin of the velocity scale, is taken to be z = 1.6408. Additional 
absorption in the CO A-X(l-O) panel at v ~ 150 kms -1 is due 
to the d-X(5-0) inter-band system. The contribution from indi- 
vidual components is shown at the position of the two reddest S i 
components (short vertical marks) by dotted and dashed profiles. 
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Observed Wavelength (A) 

Fig.3. SDSS spectrum of SDSS J 160457 +220300 fitted with the 
SDSS composite spectrum reddened by mean extinction curves 
from the Galaxy (dotted line), the LMC2 Supershell (solid line) 
and the SMC (dashed line). The arrow indicates the position of 
the 2175 A feature at z aDS = 1.64. The spectral regions used 
for x 2 minimisation are marked by solid lines at the bottom of 
the top panel. The observed spectrum is boxcar smoothed by 
10 pixels for presentation purpose. The best fit and the SDSS 
photometric points are shown in the bottom panel. 

trum and the observed QSO spectrum in the region where the 
2175 A bump is redshifted (i.e. over the observed wavelength 
range, 5203 < /1(A) < 6571). We parametrised this quantity 
by Af = (Fqso /^composite - l}- A plot of A/ versus E(B-V) to- 
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gether with the distributions of the two parameters is presented 
in Fig. [4] The positions of Q 1604+2203 in these graphs are in- 
dicated by arrows. 

As expected, QSOs with bad signal-to-noise ratio spectra 
(i.e. S/N < 10, open squares in the plots) are responsible for 
most of the scatter and a tail in the Ay distribution. Therefore, 
better confidence is achieved when selecting only high signal-to- 
noise spectra (i.e S/N > 10). Note Q 1604+2203 has S/N = 14.7. 
From the scatter plot E(B-V) versus Ay, it can be seen that the 
maximum deviation on both axis is seen for Q 1604+2203. The 
probability of getting E(B-V) > 0.15 is as low as 0.7% when we 
consider only the high S/N (> 10) spectra. This corresponds to 
a 4.7 cr significance if we approximate the E(B-V) distribution 
for spectra of S/N > 10 by a Gaussian function of mean zero and 
standard deviation 0.035. 

There is only one quasar (SDSS Jl 3 1903+43 1034) in the 
high S/N control sample with large E(B-V) in addition to 
Q 1604+2203. Interestingly, two strong Mgn systems are 
present in the spectrum of this QSO. Selecting quasar with red 
colours could therefore serve as a direct way to search for dusty 
intervening absorbers. Concerning the 2175 A absorption it- 
self, only one quasar (the one presented here) has A/ < -0.10. 
The corresponding probability is thus 0.4%. The mean A/ is 
zero with a 1 cr dispersion of 0.035. So if we assume that the 
A/ distribution for spectra with S/N > 10 is well modelled 
by a Gaussian distribution, the confidence level on the UV 
bump detection is 4.3 cr. The joint probability of measuring E(B- 
V) > 0.15 and Af < -0.10 by pure chance coincidence is 0.3%. 
In short, the exercise presented above rules out the possibility 
of the 2175 A structure being spurious at very high confidence 
level. 

It is possible to estimate N(Hi) from the reddening of 
the quasar, assuming the average LMC supershell relation be- 
tween the column d ensity of neutral gas and the extinction 
dGordon et aUl2003l) : 

MHi) „ , 

— ; — - ~ (6.97 + 0.67) x 10 21 cirT 2 , (1) 

Ay 

where A v = R v x E(B - V), and R v = 2.76+0.09 is the average 
value for the LMC2 supershell. 

We then get N(Hi) ~ 5 x 10 21 cirT 2 . We can es- 
timate a lower limit on N(Hi) using the largest value of 
A V /MH I) = 3.76 x lQ- 22 mag/(H cm" 2 ) measured in the LMC2 
bv iDobashi etaf] <|2008), JV(Hi) > 2 x 10 21 cm 2 . We note 
that using E(B-V) values obtained with different extinction laws 
(Table [2]) and the corresponding Ay /N(H i) ratios provide very 
similar results. It would be very interesting to obtain a direct 
measurement of A^(H i). However, the faintness of the quasar in 
the blue would make any determination with UVES very diffi- 
cult. 

6. Discussion 

We have reported the detection of diffuse molecular gas at z- d bs - 
1.64 towards SDSS J160457.50+220300.5, following careful se- 
lection of the target among several thousand quasars from the 
SDSS. 

High spectral resolution observations with the Ultraviolet 
and Visual Echelle Spectrograph reveal the presence of CO 
molecules in the absorbing gas. This is only the second such 
detection in an intervening DLA along a quasar line of sight. 
From modelling the CO absorption with a multiple Voigt-profile 
pattern, we measure the CO column density to be N(CO) = 




-0.2 -0.1 0.1 0.2 20 40 60 80 100 

E(B-V) Number 



Fig. 4. Distributions of reddening (E(B-V)) and strength of the 
2175 A bump (Af) for the control sample of 529 quasars 
with z em e [1.969,1.989]. The arrows mark the position 
of SDSS J160457+220300 in the different panels. Red filled 
squares and histograms indicate measurements in spectra with 
signal-to-noise ratios larger than 10. The top right panel shows 
the distribution of the signal-to-noise ratios. 

4 x 10 14 cirT 2 and an excitation temperature of the molecules, 
6 < r ex < 16 K at the 5 <x c.l., consistent with or higher than the 
expected temperature of the CMBR, Tcmbr = 7.2 K at z = 1.64. 

The observed SED of the quasar is well reproduced if dust 
reddening happens in the CO absorber with an extinction law 
corresponding to that observed in the LMC2 supershell includ- 
ing a 2175 A UV bump. We derive significant reddening (E(B- 
V)=0.27) and a depletion pattern similar to what is seen in the 
diffuse molecular medium of the local interstellar medium. 

Statistical studies have shown that DLAs as an overall popu- 
lation are not producing significant reddening of the background 
quasa r, with E(B-V) l ess tha n about 0.04 mag (Murphv & Liske, 
2004; lEllison et all l2005h . Moreover, IVladilo et all (12008) 
found a low mean extinction-to-gas ratio in DLAs Ay/N(H i)~2- 
4xl0~ 23 mag cm 2 in the SDSS-DR5. In addition, while strong 
Can absorbers seem to contain larger amounts of dust, they 
still produce on ly moderate redde ning of the background quasar 
(E(B-V) ~ O.h lWild et all 120061) . Significant reddening, simi- 
lar to that derived here, has only been found to date for a few 
individual absorpt i on sy stems at low and intermediate redshift. 
lJunkkarinen et alj d2004 report E(B-V) = 0.23 for the peculiar 
absorber at z a b s = 0.52 towards AO 0235+164 which also ex- 
hibits a 2175 A bump at the same redshift. E(B-V) = 0.16 but 
no UV bump has been found toward SD SS J1323-0021 in a 
high-metallicity sub -DLA at z abs = 0.72 dKhare et all [2004). 
lEllison et al.l 02008) detected DIBs in a Can absorber at z- d bs = 
0.16 to wards a reddened quas ar with E(B-V) = 0.23. More re- 
cently, Sriana ndet al.l d2008al) measured E(B-V) about 0.3 and 
detected 2175 A dust features in two z ~ 1.3 21-cm absorbers. 

It is interesting to note that the depletion pattern of the 
present system is similar to what has been observed in the 
z abs = 2.418 DLA towards SDSS J143912+1 1 1740 where CO 
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is also detected dSrianand et alll200 8ri). The present system has 
CO and S i column densities six times larger than in the previous 
system when integrat ed column densities of Zn n, Si n and Fe n 
are only twice larger (Noterda eme et alll2008bl) . 

The colours of Q 1604+2203 make this object laying out 
of th e quasar locus used b y SDSS to select spectroscopic tar- 
gets (Ric hards et all 120021) . More specifically, the object was 
flagged as TARGET_QSO_REJECT upon its colours and serendipi- 
tously assigned an excess spectroscopic fibre left over after the 
main samples of galax ies, LRGs, and quasars had been tiled 
dStoughton et all 120021) . We also note that this object is absent 
from the photometri c catalogue of ~ 10 6 quasars in SDSS-DR6 
dRichards et all 120091) . Similarly, the dusty absorber at z =1.3 
towards SDSS J085244+343540 dSrianand et alll2008ah . is also 
absent from this catalogue and selected as SDSS spectroscopic 
target only because the background quasar is radio-loud. This 
once again reiterates the fact that lines of sight towards colour- 
selected quasars are probably biased against the detection of the 
most relevant component of the interstellar medium, i.e., the dif- 
fuse molecular and translucent lines of sight. 

It has been four decades since the first detectio ns of the 
2175 A extinction feature in the spectra of hot stars dStecherl 
Basil. While it seems now clear that the presence of this bump is 
related to that of small dust grains, the composition of these lat- 
ter remains debated. Theoretical and experimental studies have 
shown that carbon-rich organic grains and polycyclic aromatic 
hydrocarbons (PAHs) are the most likely carriers of the 2175 A 
feature. Amorphous silicates, abundant i n interplanetary dus t 
particles could also contribute to the bump (Brad ley et all l2005). 

The absence of UV b ump in SMC lines-o f-sight and in most 
star-burst galaxies lead Gor don et al.l d 1997b to propose a trend 
in dust properties with star formation intensity. Indeed, the UV 
bump carriers could be extremely sensitive to the local chemical 
enrichment and to the energy input in the gas (e. g. lWhitteil2003t 
iGordon et all |2003HNo11 et alll2007l) . T his is supported by the- 
oretical work in the case of PAH clusters dRapacioli et all f2006) 
and could explain why such feature is ab sent in the environment 
of most gamma-ra y bursts (see however Elfasdottir et all 2009; 
iLiang & U 1200 9) and Lyman break galaxies dViih et all 2003). 
and only present when the dust is shielded enough from the in- 
cident UV radiation. 

An efficient shielding of the UV radiation field together with 
large metal abundances may thus be nece ssary for the pres - 
ence of small carbon-rich grains (see also ISloan et all [2008). 
These are the best conditions to efficiently form molecules (e.g . 
Ledou x et all 2003HSrianand et all 120051: iPetitiean et ail 2006: 
Noterdaeme et al., 2008a), and in particular CO. 

The observations reported here demonstrate that carbona- 
ceous molecules and grains already exist at z = 1 .64. This opens 
the exciting prospects of studying organic chemistry at high red- 
shift. 
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